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Abstract 

The blue-green photo-luminescence emitted by pure and electron-doped strontium titanate un- 
der intense pulsed near-ultraviolet excitation is studied experimentally, as a function of excitation 
intensity and temperature. Both emission spectra and time-resolved decays of the emission are 
measured and analyzed in the framework of simple phenomenological models. We find an inter- 
esting blue-to-green transition occurring for increasing temperatures in pure samples, which is 
instead absent in doped materials. The luminescence yield and decay rate measured as a function 
of temperature can be modeled well as standard activated behaviors. The leading electron-hole 
recombination process taking place in the initial decay is established to be second-order, or bi- 
molecular, in contrast to recent reports favoring a third-order interpretation as an Auger process. 
The temporal decay of the luminescence can be described well by a model based on two interacting 
populations of excitations, respectively identified with interacting defect-trapped (possibly forming 
excitons) and mobile charges. Finally, from the measured doping and sample dependence of the 
luminescence yield, we conclude that the radiative centers responsible for the luminescence are 
probably intrinsic structural defects other than bulk oxygen vacancies. 

PACS numbers: 
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I. INTRODUCTION 



Strontium titanate, or SrTi03 (STO), is among the most widely investigated perovskite 
oxides, owing both to its potential for novel electronic applications and to its widespread use 
as a substrate for the epitaxial growth of strongly-correlated electronic materials, such as 
superconducting cuprates or colossal-magnetoresistive manganites. Recently, the puzzling 
transport properties of its interface with other insulating oxides have also attracted much 
interest .->2 Despite this large effort, many properties of this material still await complete 
clarification. 

Intrinsic STO is a band insulator, characterized by a huge static dielectric constant re- 
sulting from the rather soft bonding of the small Ti 4+ ion to the surrounding octahedral 2 ~ 
cage. Its conduction band (CB) is composed of states having mainly Ti 3d t2 g character, 
while its valence band (VB) has dominantly O 2p character, with an upper edge located 
away from the V point in the Brillouin zone.- This results in an indirect gap of 3.2-3.3 eV, 
while the direct (optical) gap is 3.4-3.7 eV, with a non negligible sample and temperature 
dependence.-" 4 - Upon ra-type doping, generally achieved either by introducing O vacancies 
or by chemical substitution (e.g., La 3+ for Sr 2+ or Nb 5+ for Ti 4+ ), STO becomes a conduc- 
tor with a relatively large low-temperature mobility.- 1 ^ In this regime, it is known that the 
charge carriers are dressed by the interaction with the lattice and seem to behave mainly as 
large polarons, although with several nonstandard features.-^ 1 ^ When intrinsic STO is irra- 
diated with ultraviolet (UV) light, photogenerated electrons and holes (e,h) appear to both 
contribute to the material photoconductivity— Theoretical calculations would also indicate 
that holes in pure STO are not strongly coupled to phonons and keep their bare mass.— 

The photoluminescence (PL) properties of STO are possibly even more puzzling and 
controversial than its transport ones. The greenish luminescence (GL) having a maxi- 
mum at 2.2-2.4 eV of photon energy (wavelength A ~ 500 nm) that is emitted by pure 
STO at low temperature under exposure to UV or X-ray radiation has been known since 
a long time,— and is generally ascribed to the decay of intrinsic self-trapped excitons 
(STE).— A STE can be roughly depicted as a tightly bound state of a hole and a Ti 3+ 
polaron.— >2S However, this purely intrinsic scenario has been recently called into question 
by Mochizuki et al.,— >22 who argued for a crucial role of defects and possibly of surfaces in 
this GL emission. 
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Only recently, another photoluminescence emission taking place in the blue (BL), with its 
maximum at 2.8-2.9 eV (A ~ 430 nm), was reported for STO at room temperature.— 1 ^ This 
BL emission, potentially useful for optoelectronic applications, is well visible both in intrinsic 
samples at sufficiently high excitation intensities 2 ^ and in suitably n-doped samples, in the 
latter case at much lower excitation intensities.— 1 ^ A similar blue luminescence was also 
observed for intense electron-beam excitation.— 1 ^ At low temperatures, the blue emission 
is accompanied by a spectrally-narrow near-UV emission (UVL) located at 3.2 eV (A ~ 390 
nm), i.e. at band edge.— Moreover, in some cases the BL may also be accompanied 
by a long green "tail" covering a spectrum similar to the GL discussed above, but still 
visible at room temperature.— 1 ^ It is not clear what determines the appearance of this 
high temperature GL component, but surface state of oxidation seems to play an important 
role.— This GL is clearly visible at high excitation intensities.— 1 ^ 

The underlying nature of this room temperature BL is still far from clear. Although the 
BL yield is enhanced by n-doping, its spectrum looks nearly identical for different kinds 
of doping as well as for intrinsic samples under intense excitation (but this is only true 
at room temperature, as we will see below) thus pointing to a role of conduction 

band electrons rather than donor levels in the enhancement.— When studied as a function 
of temperature, in contrast with the GL which vanishes quickly above 30-50 K, the BL 
exhibits a yield that is steadily increasing with temperature, reaching a maximum at about 
160 K and then decreasing slowly (a faster decrease is however observed at temperatures 
substantially higher than room temperature, as we will see below).— ^ Studied as a function 
of the excitation ffuence (for short laser pulse excitation), the BL yield shows no saturation 
up to very high ffuences, in the mJ/cm 2 range, while the low-temperature GL has a much 
smaller saturation threshold.— 1 ^ More precisely, a detailed recent study by Yasuda et al. 
has shown that in pure STO the BL yield is actually quadratic in the excitation pulse ffuence 
up to about 1 mJ/cm 2 .— Above this value there is initially a crossover to an approximately 
linear behavior and then at 30-40 mJ/cm 2 a full saturation.— A strong dependence of 
the BL spectrum and yield on a previous surface treatment with fluorhydric acid has been 
also reported,— possibly pointing to a role of the surface or of other crystal imperfections in 
the luminescence process, although these results could be also explained as resulting from 
a variation of the surface-induced fluorescence quenching, thus unrelated with the radiative 
process itself (in Ref. |27| the excitation wavelength was 325 nm, at which the penetration 
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length in STO is of only few tens of nanometers, thus enhancing the role of the surface). 

Even more intriguing is the PL dynamical behavior studied as a function of time following 
a very short pulse excitation. It is clearly established that the PL decay does not follow a 
simple exponential behavior. At low temperatures, the GL band is associated with a slow 
power-law decay, with a very strong temperature dependence, typical of an untrapping- 
rate-limited "bimolecular" dynamics.— »i£j2I The BL dynamics has been studied in detail by 
us as a function of excitation energy (in a strong excitation regime), and we observed an 
interesting nonlinear dynamics, which can be modeled as a quite simple two-component 
decay— One component seems to be exponential, or "unimolecular" , while the other one 
can be fitted well by a bimolecular power-law decay— Although apparently similar to the 
bimolecular behavior of GL seen at low temperatures, the BL one is much faster (also at 
low temperatures) and hence must be ascribed to a different process.— The two dynamical 
components are present in the whole BL spectrum, including its green "tail" , apparently with 
no significant wavelength dependence.— Doped samples exhibit a similar two-component 
dynamics as pure ones and, at high excitation fluences, also similar decay times and yields, 
with no evident doping- induced enhancement.— However, as shown by Yasuda et al, at 
smaller excitation intensities both the overall luminescence yield and the exponential decay 
rate of the unimolecular component are found to be strongly dependent on the dopant 
concentration.— In the same paper, in contrast with our previous results, Yasuda et al. 
claim that the two-component decay is actually best described by an Auger trimolecular 
process acting together with the unimolecular one, thus leaving the issue of the leading 
recombination process governing the PL decay undecided. Properly assessing the strength 
of the Auger recombination may be important for proposed applications of STO in opto- 
thermionic refrigeration ."22, We will come back to this issue below. 

In this work, we investigate further the physics underlying the blue luminescence of STO 
by analyzing the PL spectra and temporal decays, for both intrinsic and doped samples, as a 
function of excitation energy and temperature. The temperatures of our measurements are 
high enough to make the activated luminescence quenching evident. Moreover, we compare 
in detail the predictions of different models for the PL dynamical behavior, in order to 
establish the nature of the recombination processes involved in the decay and to shed light 
on the underlying electronic mechanisms. 

This paper is structured as follows. Section 2 describes the experimental procedures. The 



4 



wavelength-resolved studies as a function of excitation energy and temperature are reported 
and analyzed, within a simple model, in Sec. 3. The corresponding time-resolved studies are 
then reported in Sec. 4, together with a first phenomenological modeling. In Sec. 5 we tackle 
the question of a more detailed modeling of the dynamical behavior. The final Section 6 
includes a discussion of the physics underlying the dynamical models used for interpreting 
the PL decay and a summary of our main results. 

II. EXPERIMENTS 

The samples used in this work were of the following three kinds: (i) five stoichiometric 
intrinsic (lOO)-oriented STO single crystals (I-STO), 1x5x5 mm 3 in size, produced by 
four different companies (SurfaceNet GmbH, CrysTec GmbH, Crystal GmbH, eSCeTe B.V.) 
by the flame-fusion Verneuil method, with specified impurity levels all below 150 ppm, 
and used as received; (ii) two samples of Nb-doped STO crystals (N-STO), with a Nb molar 
concentration of 0.2%, and having the same orientation and geometry as the I-STO samples; 
(iii) one sample doped with oxygen vacancies (O-STO), obtained from a pure sample by 
annealing for 1 h at 950°C and 10 -9 mbar (base pressure 10 -11 mbar). While I-STO samples 
look transparent and are verified to be good insulators, both N-STO and O-STO samples 
look black/dark-blue opaque and are conducting. In O-STO, from resistivity measurements 
performed on a thin film annealed by the same procedure, we estimate an induced carrier 
density of 3 x 10 17 cm -3 , or about 0.002%, rather small but still much higher than typical 
residual oxygen vacancy concentrations in nominally stoichiometric samples.— I-STO and 
N-STO samples did not show aging or hysteretic behavior due to oxygen exchange with 
atmosphere, as we checked by repeated measurements. Instead, O-STO samples turned 
transparent again when heated above 250 °C in air, clearly showing re-oxidation of the 
vacancies. For this reason, all our temperature behavior studies were limited to I-STO 
and N-STO samples. During measurement, the samples were held into a thermostat for 
temperature control to within 0.1 K (with an accuracy within few K). The temperature was 
scanned between 300 K and 900 K. 

In all our experiments, the excitation was induced by 3.49 eV UV photons (A = 355 nm) 
in 25-ps-long laser pulses at 10 Hz repetition rate focused to a gaussian spot having a radius 
of 1.2 ± 0.1 mm at 1/e 2 of maximum. The energy per pulse was varied from 40 fiJ to 2 mJ, 
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corresponding to an excitation fluence U ranging from a minimum of 2 to a maximum of 100 
mJ/cm 2 (value at the spot center, corresponding to twice the spatial-average value). This 
goes up to much higher values than those investigated by others, including Mochizuki et al. 
(up to 14 mJ/cm 2 )2i and Yasuda et al. (up to 10 mJ/cm 2 } 3 -^. For our highest fluence of 100 
mJ/cm 2 , taking into account the 25% reflection and assuming an optical penetration length 
of about 1 /im, 4 we estimate a peak density of photogenerated e-h pairs as high as 1.2 x 10 21 
cm -3 . This corresponds to a fluence-to-density conversion factor (FDCF) a = 1.2 x 10 22 
cm _1 /J. It must be noted, however, that the UV penetration length is highly uncertain, as 
different samples have shown fairly different absorption edges in previous reports. We stress 
that, despite our large excitation fluences, no visible photoinduced damage of the sample 
surfaces was induced during our experiments and we never observed irreversible variations 
of the signal ClS db function of excitation intensity. 

The luminescence emitted from the sample was collected by a lens system imaging the 
illuminated sample spot onto the detector head, after blocking the (much stronger) elastic 
scattering by a long-pass filter with a cutoff wavelength of 375 nm. When recording the 
luminescence spectra, for detection we used a grating-monochromator and a photomultiplier 
and integrated the signal in time (with a 50 ns time gate). In time-resolved measurements, 
the luminescence was instead detected with a photodiode (PD) having a rise-time of about 
150 ps. In time-resolved experiments the entire luminescence spectrum was integrated. The 
PD signal was acquired by a 20 Gsample/s digital oscilloscope having an analog bandwidth 
of 5 GHz. The response function r(t) of this apparatus was acquired by measuring the signal 
given by the elastic scattering of the excitation pulse (taken after removing the long-pass 
filter). 

III. MEASUREMENT RESULTS: SPECTRA 

The data shown in Fig. [1] refer to a typical I-STO (panels a and b) and a typical N-STO 
sample (panels c and d), respectively under 2.2 (panels a and c) and 22 mJ/cm 2 (panels 
b and d) excitation fluences, respectively. It is seen that the overall shape of all spectra 
is asymmetrical, with both GL and BL emissions present and clearly distinguishable at 
all temperatures, in spite of a substantial overlap. At a closer inspection, an additional 
spectral contribution apparently located at about 3.0 eV (A ~ 400 nm, in the violet), 
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FIG. 1: (Color online) Photoluminescence spectra from I-STO (panels a and b) and N-STO (panels 
c and d) samples at various temperatures and at an excitation fluence of 2.2 mJ/cm 2 (panels a and 
c) or 22 mJ/cm 2 (panels b and d). Labels GL, BL, and VL (green, blue, and violet luminescence) 
in panel (a) mark the peaks of the three separate luminescence bands that can be singled out in 
our spectra. A blue-to-green spectral transition (a red-shift of the peak) of the photoluminescence 
is seen in the I-STO spectra, while it is absent in the N-STO doped samples. The lines in panel 
(b) are best-fit curves based on our theory. The lines in the other panels are guides to the eye. 

henceforth called VL (see panel a), can be singled out, particularly in I-STO samples (it is 
however possible that this emission is actually the same as the UVL band mentioned in the 
Introduction, after a reshaping due to the detection-line long-pass filter used in our setup). 
This VL appears as a small shoulder in the lower temperature measurements, while it is more 
clearly separated, although strongly depressed, at higher temperatures. Our spectra are 
qualitatively consistent with literature data. Three distinct emissions with similar spectral 
position were reported in Ref. 



STO reported in Ref. 
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though at lower temperature. The spectrum of La-doped 
closely resembles that of our N-STO samples at 300 K. However, 
compared to previous data, ours seem to show a comparatively higher yield in the GL region. 
In all samples the GL contribution appears to increase more weakly than the BL one for an 
increasing excitation intensity, as can be seen in Fig. CD by comparing panels (a) with (b) 
and (c) with (d). Therefore, the GL emission is probably more saturated than the BL one, 
pointing to the presence of different classes of emitters as responsible for the two bands. 

As a function of temperature T, the most striking feature seen in the PL spectra is a strong 
red-shift of the BL maximum occurring for increasing T in I-STO samples. I-STO samples 
at high temperature emit almost exclusively in the green, with only a small residual emission 
in the blue (which we ascribe to the VL band) . This blue-to-green thermal transition of the 
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luminescence can be clearly seen even with the naked eye. The fact that this shift does not 
occur (or is very small) in N-STO doped samples, makes the spectra of I-STO and N-STO 
samples become very different at high temperatures. This is an important observation, as 
it is the first clear qualitative difference seen in the luminescence spectra of doped samples 
as compared to pure ones. In contrast to what stated before,— 1 ^ this shows that the role of 
doping is not limited to introducing additional charge carriers in the system, but somehow 
affects also the properties of the radiative and/or non-radiative recombination centers. 

In order to analyze the data quantitatively, we will refer to a specific model for the 
spectral shape I(ui) of each component of the luminescence spectrum (GL, BL, and VL), 
treating them as independent. As previously stated, the low-temperature GL was interpreted 
in a quantitative way by Leonelli and Brebner in terms of the annihilation of self-trapped 
excitons of given energy Eq.— In this picture, the broadening of the band is due to the 
random emission of several optical phonons, each of energy hfl ~ 90 meV, giving rise to the 
following expression: 

I(u>) ~ I e- s ° V Sgr/(27rn!) 

whose key parameter is the Huang- Rhys factor Sq, related to the strength of the electron- 
lattice interaction and fixing the average number of emitted phonons per recombination. T 
is the width of the band associated with a given n-phonon process, taken to be Lorentzian. 
Although devised for self-trapped excitons, Eq. §T§ applies equally well to the case of defect- 
assisted recombinations - in which defects mediate the coupling of the electronic excitation 
to the lattice - , so it may be actually considered as a semi-phenomenological model that 
can describe different microscopic scenarios. Leonelli and Brebner found that a value as 
high as So = 5.7 produced a nice fit to the data, meaning that the maximum of intensity at 
2.4 eV is well shifted with respect to the intrinsic exciton energy, taken at Eq ~ 2.9 eV. We 
plotted the I(oS) resulting from Eq. (JTJ, keeping the same quoted values of the parameters Sq 
and Eq and assuming a broadening factor T = 0.12 eV: even with no adjustable parameters 
(except for the overall amplitude scale), we could describe in this way quite accurately the 
GL band shape in the region where it is well separated from the other bands (see, e.g., the 
highest temperature spectrum in Fig. [J}}). Hence, we assumed that the GL band keeps the 
same form in all spectra, except for a temperature dependent scale factor, and turned to the 
problem of the BL and VL bands. Tentatively, given its general applicability, we used Eq. 

8 




FIG. 2: (Color online) Temperature behavior of the amplitudes of the three spectral bands GL 
(triangles), BL (circles), and VL (squares) as seen in the I-STO (panel a) and N-STO (panel b) 
photoluminescence. The excitation fluence was 2.2 mJ/cm 2 . Solid lines are best fits based on 
Eq. §2$; the resulting activation energies are reported in the legend. The case of high excitation 
fluence (22 mJ/cm 2 ) gives similar results, except for a slightly smaller activation energy of the VL 
component (« 0.2 eV instead of ~ 0.3 eV). 

([I]) also for these bands, and after suitable adjusting of the parameters E and So we could 
achieve a satisfactory fit of all data. The phonon energy hQ and the irrelevant broadening 
factor T were always kept fixed at the values 88 meV and 0.12 meV, respectively, without any 
attempt of optimization. The optimal values of the characteristic energies and of the Huang- 
Rhys factors for the GL, BL and VL components were determined by a global fit procedure 
as ~2.9 eV, ~2.9 eV, ~3.0 eV, and ~5.5, ~1, ~1, respectively, independently of fluence and 
temperature, and allowing for only a slight sample dependence. In addition, we used as free 
fit parameters the respective amplitudes Iq of the three bands at each temperature. Typical 
results for pure and doped samples are reported in Fig. [2j As a last step, we observed 
that allowing for a slight decrease of the characteristic energy of the VL contribution for 
increasing temperature (in any case below 0.1 eV) the fit was further improved. This shift 
can be probably associated with the known temperature dependence of the STO gap.— An 
example of the fit results is given in Fig. [T]d. The overall fit quality is as good in all I-STO 
and N-STO investigated samples. 

Figure [2^ shows that the blue-to-green transition of I-STO samples can be explained as a 
thermal quenching of the VL and BL components occurring at lower temperatures than for 
the GL one. This is not the case of N-STO samples, for which the three amplitudes have a 
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FIG. 3: (Color online) Histogram of the PL yield (panel a) and PL-tail exponential decay time r e 
(panel b), for different samples at room temperature. Samples 1-5 are I-STO (blue online), sample 
6 is O-STO (red online) and samples 7-8 are N-STO (green online). 

more similar behavior with temperature, as shown in Fig. 0x The thermal quenching of the 
PL amplitudes Iq(T) can be approximately modeled by the following standard Arrhenius 
activation law: 



J„(T) oc r(T) oc \a + be~ Ea/{k » T) 



(2) 



where r is the characteristic PL decay time, E a an activation energy giving the potential 
barrier of the competing non-radiative relaxation channels, kf, the Boltzmann constant and 
a, b are constants that give the relative weight of the temperature independent (typically 
radiative) and thermally activated (typically non-radiative) contributions to the decay, re- 
spectively. The best fit results based on Eq. (j2J) are also shown in Fig. [2j The corresponding 
activation energies are all of the order of tenths of eV, with somewhat smaller values for 
GL (0.16-0.18 eV), and slightly larger for BL and VL (0.23-0.29 eV), with the exception of 
the BL in I-STO, which is found to be about 0.8 eV. According to these best-fit results, the 
stronger thermal quenching of the BL and VL bands compared to GL is not associated to 
a smaller activation energy (it is actually larger), but to a relatively larger weight of the 
non-radiative channels with respect to the radiative ones. 

Finally, the spectrally-integrated PL yield at a pump fluence of 2.2 mJ/cm 2 as a function 
of the sample at room temperature is shown in Fig. [3^,. The most important things to notice 
here are the following: (i) there is a significant sample-to-sample dependence of the yield 
in nominally identical I-STO samples; (ii) doped O-STO and N-STO samples do not show 
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FIG. 4: (Color online) Time decay of the PL signal of a I-STO sample following its excitation by 
a UV picosecond pulse, for different excitation fluences, at room temperature. Data are shown 
as gray dots (blue online). The solid lines are the result of a global best fit based on our models 
(including a final convolution with the measured instrumental response time r(t)). The black line 
is based on the C2PUBvl bimolecular model and the gray line (red online) is based on the 1PUT 
trimolecular model. Data and curves referring to different fluence values are vertically shifted for 
clarity. 

any enhanced yield, despite the presence of many additional donor defects. We will discuss 
these findings in Sec. 6. 



IV. MEASUREMENTS RESULTS: TEMPORAL DECAYS 



A typical set of PL temporal decays measured from a I-STO sample for various excitation 
pulse fluences U is shown in Figs. HH5J The case of N-STO is very similar (an example is 



reported in Fig. 2 of Ref. |29| ). As already noted in our previous works and confirmed by 
Yasuda et aZ.,— -M-^L the PL initial decay becomes faster for higher excitation energies, while 
the final part of the decay (the "tail") varies only in its amplitude, relative to the peak, 
but not in its rate. One can therefore phenomenologically single-out two distinct regimes: 
the initial fast decay, with an excitation-dependent characteristic decay time, and the final 
tail, with a well defined, excitation-independent, exponential decay time r e . In I-STO, the 
fast decay is typically in the range 1-2 ns (for our excitation fluences) and is only weakly 
sample dependent, while the slow tail time constant ranges from 6 to 23 ns, depending on 
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FIG. 5: (Color online) Time decays of the PL signal of a I-STO sample at room temperature, for 
different excitation fluences, shown here in semi-logarithmic scale to highlight the exponential PL 
tail. Data are shown as dots, while the solid lines are the result of a global best fit based on the 
C2PUBvl model. The wiggles seen in the tail for both data and model predictions are due to the 
detection-system response function r(t), which is shown in the inset (inset axes labels and units 
are the same as for main panel). 

the sample, as shown in Fig. [3h>. In N-STO the fast decay rate is of the same order as in 
I-STO, while the exponential decay becomes somewhat faster (r e w 5 ns), as an effect of 
doping. 31 

In Figs. HJl5] the decay signals are normalized to their maximum, for clarity. However, 
also the PL signal amplitude varies strongly with the excitation intensity. More precisely, 
as mentioned in the Introduction, we observed an approximately linear dependence of the 
time-integrated PL signal, i.e. of the overall PL yield, on excitation fluenc e up to 30-40 
mJ/cm 2 , while for even higher fluences we find a saturation (see Fig. 2 of Ref. 



of Ref. 



23 and Fig. 3 



29 ). The PL signal maximum typically shows a mixed linear-quadratic behavior 



with excitation fluence, with a relatively large scatter of the data, for reasons which we have 
not yet identified. The exact threshold for saturation is also sample dependent, pointing 
again to an important role of the defects in the PL radiative channels. Moreover, the linear 
yield behavior we find in our data is actually already a partly saturated one, as for lower 
excitation fluences a quadratic behavior is observed instead.— 

To specify a characteristic experimental decay time for each given PL signal, we use 
its full- width- at- 1/e-of-the- maximum (FWl/eM) time ty e . Equivalently, its inverse ty e 
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FIG. 6: (Color online) Experimental decay rates as a function of excitation fluence and temperature, 
for a I-STO sample. The decay rate is here defined as the inverse of the FWl/eM decay time 
t\j e . The upper x axis gives our estimate of the density of photoinduced electron- hole pairs after 
excitation. The lines are linear best fits of the fluence dependence. 

provides a characteristic decay rate. The measured decay rates ty e versus excitation fluence 
U, for different temperatures, in a I-STO sample, are shown in Fig. [6j N-STO samples give 
very similar results. In the investigated range, this decay rate shows an approximate linear 
dependence on the excitation fluence, as shown in the figure. However, since the measured 
decay times ty e are close to the characteristic response time of our detection apparatus, these 
raw data for tu e are significantly larger than the actual decay times of the luminescence. 
To take care of this, in principle we should deconvolve the measured decay signal and the 
response function of our setup. However, the numerical deconvolution of noisy data is 
known to be problematic, so that, following the standard procedure, we used the inverse 
approach. Given a theoretical model I(t) for the PL decay containing the characteristic 
times to be determined as adjustable parameters, we first convolve it with the measured 
response function r(t), thus obtaining a predicted signal S(t) = (r * I)(t). The latter is 
then compared with the measured signal, thus finding the best- fit values of the adjustable 
parameters. In this way, the best-fit values of the characteristic decay times appearing as 
parameters in the model function I(t) will correspond to actual PL decay times, without 
significant distortions due to the instrumental response function. The final time-resolution 
that can be achieved with this approach is limited only by the signal-to- noise ratio of our 
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FIG. 7: (Color online) Characteristic PL decay times t\ = r e (squares) and 7~2 (circles), the latter 
at the reference fluence of 10 mJ/cm 2 , versus temperature T, for a I-STO sample (panel a) and a 
N-STO sample (panel b). The lines are best fits based on the activated behavior given by Eq. ([2]). 
The resulting activation energies are given in the legend. 

data (typically 10 2 — 10 3 ). In this Section, we wish to analyze the data without relying 
on too specific physical assumptions: ideally, we would like to adopt a model-independent 
description, postponing the analysis of specific models to the following Section. Therefore, 
we adopt here the "phenomenological" model already used in Ref. |28| (corresponding to the 



2PUB considered in next Section). According to this model, the PL intensity as a function of 
time is given by the sum of a pure unimolecular exponential decay with time constant T\ and 
a pure bimolecular decay with excitation-dependent time constant r 2 = 1/(70:2?/), where 7 
and «2 are constants.— The fast FWl/eM time t\/ e will be approximately proportional to 
the value of r 2 . For each given sample and temperature, we may assign a typical value to the 
time constant r 2 by fixing a reference excitation fluence U\ at which it must be evaluated. 
We choose here a value for U\ that is in the order of our experimental range, i.e. U\ = 10 
m J / cm 2 . 

The best-fit values of the exponential slower decay time T\ = r e and of the faster initial 
decay time r 2 (at the reference fluence of 10 mJ/cm 2 ) for a I-STO and a N-STO samples 
are shown in Fig. [7] (the fit procedure is described in more detail in the next Section). As in 
the case of the PL spectral amplitudes, all decay time temperature behaviors can be fitted 
fairly well by the Arrhenius law given by Eq. ([2]). The resulting activation energies are given 
in the figure legend and are in the range 0.2 — 0.3 eV, of the same order as those obtained 
from the behavior of the PL band amplitudes. 
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FIG. 8: (Color online) Photoluminescence yield, computed from the time-integral of the decay 
signal, as a function of temperature, for a I-STO and a N-STO samples. Solid lines are best fits 
based on the activated behavior given by Eq. ([2]). The resulting activation energies are given in 
the legend. 

To conclude this Section, we report in Fig.[8]the measured temperature behavior of the PL 
yields obtained from the time integral of the decays, together with the usual activation- law 
best-fits. The resulting activation energies are again roughly consistent with the dynamical 
ones and with the spectral amplitude ones. No significant dependence on doping is found 
in this case. 

V. MODELING THE TEMPORAL DECAY 

Although it is firmly established that the initial decay rate of the PL at sufficiently 
high excitation intensities behaves nonlinearly,— i2iL2i the exact decay law and the underly- 
ing microscopic relaxation mechanisms are currently controversial, as we mentioned in the 
Introduction. In this Section we introduce and compare several different models for the PL 
decay, listed in Table I, aiming at identifying the most effective one in describing our data, 
which can in turn offer indications about the underlying microscopic physics (to be discussed 
in the next Section). In particular, we are especially interested in assessing the order of the 
leading recombination mechanism, e.g., bimolecular versus trimolecular. 

Each tested model is labeled with a code (first column in Table I) which synthesizes the 
most important assumptions on which it is based: (i) the figure before the letter "P" gives 
the number of decaying electron populations considered in the model (either 1 or 2), with 
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a "C" denoting the more specific case of coupled populations; (ii) the letters "U, B, T" 
stand respectively for unimolecular, bimolecular and trimolecular recombination mechanism 
(a unimolecular recombination process is taken to be present in all models); "vl" and "v2" 
denote different variants of the same basic model. 

In this work, we adopted a "global fit" procedure, that is we fitted simultaneously all the 
decays measured for a given sample at a given temperature but for different values of the 
fluence U, for a single choice of the adjustable parameters. For each model and for all samples 
and temperatures, we have performed global best fits of the following two kinds: (i) "without 
the amplitudes" , i.e. on signals which had been previously normalized to their maximum; (ii) 
"with the amplitudes" . The reason for using both methods is that the PL decay amplitudes 
were usually subject to a significant scatter, not present in the decay functional form, and for 
the highest excitation fluences also to some degree of saturation, which may not be properly 
taken into account in our simple models. Therefore, approach (i) is more appropriate in 
order to focus on the capability of our models to predict the ^-dependence of the PL decay 
functional form, and in particular of the PL decay rates, regardless of the amplitude behavior. 
Approach (ii) tests the models in their predictive power for both PL rates and amplitudes, 
but tends to weight more the amplitude behavior, as it gives the strongest data variations. 
All compared models have four adjustable parameters in approach (ii), while in approach 
(i) all models have three parameters except for one (C2PUBvl), which has four. For each 
approach, we quantified the performance of each model with the fit \ 2 (normalized to the 
data variance, estimated using the measured noise before the excitation pulse), averaged 
over different samples and repeated measurements as specified in Table I. 

The first two models in Table I (1PUB and 1PUT) are based on the assumption that a 
single-population density N(t) of decaying electrons and holes (balanced in number) suffices 



for capturing the recombination dynamics revealed by the PL decay. Many past ana 



yses of 



picosecond-laser-induced PL temporal decays in semiconductors (see, e.g., Refs. I34ll35lj 36l). 
as well as the recent work of Yasuda et al. on the STO,— interpret the PL decay by a similar 
single population model. While the unimolecular and bimolecular terms are fairly common 
in solid state luminescence, the appearance of a third-order trimolecular term, typically 
ascribed to Auger processes involving two electrons and a hole or two holes and an electron, 
is only seen at very high excitation densities in indirect band-gap materials. In principle, it 
is quite reasonable to expect this behavior in the case of STO, and indeed the 1PUT model 
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TABLE I: Different PL decay models tested in this work. The two best-fit values of x 2 reported 
in the last two columns are averaged over different samples and repeated measurements, after nor- 
malizing to the minimum x 2 value obtained among all models for each given sample/measurement 
(in order to weight all samples equally). The first x 2 ( n o ampl.) is computed for decays normalized 
to their maximum, so that the behavior of the decay amplitude with excitation fluence does not 
enter the fit and the model testing is focused on the decay rates. The second x 2 (with ampl.) 
is instead computed taking also the decay amplitudes into account. The reported differences are 
statistically highly significant (the formal likelihood ratio between the best model and the others 
is of thousands of orders of magnitude). 



model code 


rate equations 


initial conditions 


radiative term 


X 2 no ampl. 


X 2 with ampl. 


1PUB 


dN/dt = -N/ti - 77V 2 


N(0) = aU 


/ = Q-yN 2 


2.1 


1.9 


1PUT 


dN/dt = -N/n - C A N 3 


same as above 


I = Q-yN 2 


3.2 


1.8 


2PUB 


dNi/dt = -N 1 /t 1 
dN 2 /dt = -7/Vf 


iVi(O) = aiU 
Ar 2 (0) = a 2 U 


I = Q 1 N 1 /r 1 +Q 2l N 2 


2.2 


1.7 


2PUT 


dN-i/dt = -N 1 /t 1 
dN 2 /dt = -C A N% 


same as above 


I = Q 1 N 1 /r 1 +Q 2 yN 2 


3.1 


1.7 


C2PUBvl 


dNi/dt = -JVi/ri - 71^1^2 
dN 2 /dt = -72 N 2 


same as above 


I = Q \dNi/dt\ 


1.0 


1.7 


C2PUBv2 


same as above 


^1(0) = a + aiU 
jV 2 (0) = a 2 U 


I = Q lr N 1 N 2 


1.3 


1.1 



is that assumed in Ref. 3l|. However, we find that the 1PUT model is very ineffective in 
describing our data. An example of the unsatisfactory results of a global fit based on the 
1PUT model is shown in Fig. H] (gray lines, red online), while the average best-fit x 2 values 
relative to the best-performing models are reported in Table I. 

We stress that by using a global (simultaneous) fit on all the decays obtained at different 
excitation energies U, we have put the models to the test not only on their capability to 
predict the single-decay functional form, but also on their capability to predict the overall 
dependence of this functional form on the excitation fluence U. This dependence is very 
sensitive to the order of the highest nonlinear term in the decay equations and therefore 
matching this dependence provides a much stricter test than matching the single-decay 
behavior. More precisely, the initial rate of decay taking place shortly after excitation is 
predicted to depend much more strongly on the excitation fluence U for the 1PUT model 
than for a second-order model, such as 1PUB (see Table I). Indeed, assuming that the PL 
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FIG. 9: (Color online) Inverse decay times versus excitation fluence for the four PL decay curves 



reported in Fig. 1 of Ref. 
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j. The times have been estimated graphically, from the initial decay 
slope in the semi-logarithmic chart. The line (red online) is a linear best-fit. The good linearity of 
the decay rates vs fluence confirms a bimolecular behavior for the initial fast decay. 

signal is generated by a given power p of the population density, I{t) oc N p (t) (e.g., p = 2 
for a bimolecular radiative recombination), the initial logarithmic decay rate for a single 
population model including both bimolecular and trimolecular terms is given by 



d\nl 



dt 



p dN 



t=o+ 



N dt 



t=o+ 



p -+ 1 N(0 + ) + C A N 2 (0~ 

Ti 

p(—+ jail + C A a 2 U 2 

VTi 



(3) 



Thus, for the 1PUT model one would expect a quadratic dependence of this initial decay 
rate on excitation fluence U, while for a second-order (bimolecular) model, having C A ~ 0, 
one expects a linear dependence. This initial decay rate is well estimated by the inverse 
FWl/eM time £7/ e . We indeed find a linear behavior, as shown in Fig. [61 and this supports 
a bimolecular model. To investigate also the hypothesis of a crossover from a trimolecular 
behavior taking place for the lower fluence range studied by Yasuda et al. (0 — 10 mJ/cm 2 ) 



to a bimolecular one in our 
extracted from Fig. 1 of Ref. 



ligher range of fluences (10 — 50 mJ/cm 2 ), we graphically 
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the initial logarithmic slope of the four reported decays. 
These values are plotted in Fig. [9] versus the excitation fluence, showing that the behavior is 
again perfectly linear in U as expected from a bimolecular model (not necessarily a single- 
population model), and hence not consistent with a third-order model such as 1PUT. We 
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note that the fluence dependence of the measured decay rates can be altered by the finite 
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were taken with a streak camera, 



response function of the apparatus. The data from Ref. 
which typically has a very fast response in the picosecond range, so that they should not 
be affected by this problem. Our data are instead affected by the slower response time of 
our equipment. Nonetheless, even after convolution with the response function, the 1PUT 
third-order model predicts a stronger dependence of the initial decay rate on the excitation 
fluence than what seen in our data, as shown already in Fig. HI and more explicitly in Fig. []JJ 

In synthesis, we can conclude that the initial faster decay is a second-order, or bimolecu- 
lar, recombination process and no significant trimolecular effect is detected in our data. We 
note however that even at fluences much smaller than ours, indirect semiconductors usually 
exhibit a strong Auger-like third-order decay.— Presumably, third-order Auger interac- 
tions are depressed in STO by its relatively large band-gap 37 and high dielectric constant, 
while bimolecular recombinations might be favored by its comparatively large density of 
intragap trapping states and by its stronger electron-phonon interactions. 

Although it behaves much better than the trimolecular model, we see from Table I that 
also the bimolecular single-population model 1PUB is not fully satisfactory in describing our 
data. For this reason, we decided to consider models based on two dynamical populations, 
N±(t) and AT 2 (t), representing for example free and trapped charges. The simplest model 
of this kind is the one already considered in our previous papers,— labeled as 2PUB, which 
corresponds to the case of two independent populations, one decaying with a unimolecular 
process and the other with a bimolecular one (see Table I). The best-fits obtained using this 
model are fairly good, but again not fully satisfactory. Moreover this model does not lend 
itself to a simple and plausible physical interpretation. However, this 2PUB model has the 
advantage of having two well separated terms describing the initial faster decay and the final 
slower tail. Therefore, it is particularly apt to describing phenomenologically the data, for 
characterizing the faster and slower decay rates of our data in a roughly model-independent 
way, as we have done in Sec. 4. For the sake of completeness, we also considered a two- 
population unimolecular + trimolecular model (2PUT), which however can be discarded 
after comparison with data (see, e.g., Fig. [TO]) . 

To go beyond the 2PUB model while remaining in the framework of a two-population 
model, it is therefore necessary to assume some kind of coupling between the two populations 
(C2PUB models). The simplest choice is to include a term proportional to the cross product 
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FIG. 10: (Color online) Inverse decay rates versus fluence as measured for a I-STO sample (squares) 
and as predicted by different models (lines), after a global best fit on the measured decays. Besides 
ours, the Yasuda et al. data already given in Fig. [UJ) are also plotted (circles). The two sets of data 
do not merge into a single smooth behavior because of the much slower response of our apparatus 
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The predictions of the C2PUBvl model 



as compared to that used by Yasuda et al. in Ref. 
(black lines), after a global fit on our data only, explain both our data (solid line, obtained after 
a convolution with our instrumental response function) and the Yasuda et al. data (dot-dash line, 
obtained from the C2PUBvl model assuming an instantaneous instrumental response). Dashed 
and dotted curves correspond to the predictions of the 1PUT model (gray line, red online) and 
of the 2PUT model (light gray line, orange online), after convolution with our response function 
(dashed lines) or for an instantaneous response (dotted lines). The gray solid line (cyano online) 
is a linear best-fit to our data. 

N\N 2 in one of the rate equations, i.e. a recombination process of one population that is 
stimulated (or assisted) by the other population. This may arise from a variety of processes, 
as will be discussed in the next Section. The solution to C2PUB rate equations (with the 
more general initial conditions given in the table for C2PUBv2, see below) is the following: 

(a + aiU)e' t/Tl 



N x {t) 
N 2 (t) 



(1 + t/T 2 p 1 ^ 2 

a 2 U 



(4) 



1 + t/r 2 

with r 2 = l/{'~y 2 <y 2 U). The model must be now completed with an assumption about the 
radiative terms. We consider here two different possible choices for this assumption, leading 
to the two model variants given in Table I. The first (C2PUBvl) is obtained from the 
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assumption that both recombination terms in the N\ population are partly radiative, with 
the same quantum efficiency Q. With this assumption, the predicted PL decay is given by 
the following expression (in which we have also set a = 0, as otherwise we would get an 
absurd nonzero PL for zero excitation): 



I(t) = Q 



e -t/Tl 



(1 + t/r 2 )^ 2 
e t/Tl 

(i + t7^) 1+71/72 



dt 

+k2U 2 -— 7 - T —- Ci) 



where k\ and k 2 are constant amplitudes, and the following relationship holds between the 
parameters: 

71 _ k 2 r 2 U _ k 2 

72 fan &iTi7ia 2 
From our best fits, we obtain a typical ratio 71/72 = 0.1 — 0.3. 

This model is in very good agreement with our data (for normalized decays), as shown 

for example in Fig. H] and [Till Fig. [10] also shows that the C2PUBvl model, after fixing 

its parameters to those giving a best fit to our data, predicts well without any further 

adjustment also the decay rates measured by Yasuda et al. for a smaller excitation fluence 

range. 31 The advantage of the C2PUBvl model, compared to the previously mentioned ones, 

is also quantitatively reflected in the normalized x 2 °f the fit, which is smaller than for the 

other models by a factor two/three, statistically very significant. However, when comparing 

normalized decays, this model has one adjustable parameter more than the others (four 

against three), so that there is still margin for doubts. Moreover, when fitting also the PL 

amplitudes this model does not perform much better than the others. This is probably 

revealing some saturation behavior of the amplitude data that is not well captured by this 

model. 

Let us now consider the second variant of the C2PUB model (C2PUBv2) that is obtained 
by the following assumption: the radiative emission is now taken to arise only from the cou- 
pling term NiN 2 , i.e. from the decays involving both populations 1 and 2. Mathematically, 
this is equivalent to setting fa = in Eq. (jSJ). The number of adjustable parameters is 
thus back to three in the case of normalized decays. Nevertheless, this C2PUBv2 model is 
almost as effective as C2PUBvl in the global best fits on normalized decays. In the best-fits 
including also the decay amplitudes (i.e. in approach ii), as can be seen from the \ 2 values 
given in Table I, this model is actually giving by far the best results (using the same num- 
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ber of adjustable parameters as for other models) provided that we use the modified initial 
conditions indicated in Table I with a ^ 0. This corresponds to a saturated dependence 
of iVi(O) versus U, and mathematically leads to the replacement k 2 U — > {k%-\U + k 2 2U 2 ) 
in Eq. (jSJ) [since we know that 71 <C 72, we also set 71 w in Eq. (jSj), in order to keep 
the number of adjustable parameters to four]. This modification of the initial conditions 
does not affect the best-fits on normalized decays because it alters only the predicted PL 
amplitudes but not the decay rates and functional forms. As we will discuss in the following 
Section, this modified choice of initial conditions lends also itself to a simple and plausible 
physical interpretation. 

Before concluding this Section, we note that for T\ 3> t 2 , both models C2PUB (vl and v2) 
predict an initial time decay of the PL of the form I(t) ~ 1 / (1+t / V 2 ) ( ' 1+71 which, for small 
values of the ratio 71/72 such as those found in our best fits, is very similar to that predicted 
by trimolecular models, i.e. I(t) ~ 1/(1 + £/t 2 ). The other bimolecular models 1PUB and 
2PUB that we analyzed predict instead an initial behavior that is inversely quadratic rather 
than linear in the time, i.e. I(t) ~ 1/(1 +t/r2) 2 . This fact probably explains why Yasuda et 
al, when fitting the single decay curve, found a better fit for their decays using a trimolecular 
model rather than a single population bimolecular one. On the other hand, our global fits 
discriminate much more effectively among the various possible models. 

VI. DISCUSSION AND CONCLUSIONS 

We now turn to discussing the microscopic physics that may underly the STO blue 
photoluminescence phenomena we have described. The photo-generated carriers may in 
principle populate different excited states, which can be both localized (trapped) or extended 
(mobile). Mobile charges will belong to the STO CB as electrons and to the VB as holes, 
in both cases probably with some degree of phonon dressing (large polarons). Localized 
charges may in principle be self-trapped and fully intrinsic (small polarons, or self-trapped 
excitons if paired) or associated with intrinsic crystalline disorder or defects such as oxygen 
vacancies, dislocations, or possibly surface states. The indirect nature of the STO bandgap 
forbids CB-VB direct recombinations, so that electron coupling to phonons, presumably 
enhanced by defect- or self-trapping, is an essential element of the luminescence process. 

Let us start by discussing the spectral and yield features of the PL, in order to identify 
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the nature of its radiative centers and of the competing non-radiative channels that may 
contribute to its quenching at high temperature. According to our fits of the PL spectra, the 
typical total energy released by the annihilation of an e-h pair is ~2.9 eV, partly dissipated 
in phonons and partly radiated. The red shift of the GL spectral component with respect to 
the BL and VL components is attributed to the larger fraction of vibrational energy released 
in the former case, implying a stronger electron-phonon coupling. The radiative processes 
giving rise to the PL can be in principle fully intrinsic, i.e. characteristic of a perfect STO 
crystal, such as phonon-assisted CB-VB direct recombinations or STE annihilation, or again 
associated with intrinsic structural lattice defects such as those mentioned above.— Extrinsic 
defects such as chemical impurities are likely to be excluded, instead, because the PL yield 
does not present the very large sample-to-sample fluctuations that are typical of impurity- 
associated luminescence (see Fig. |3]).— Moreover, the lack of any significant doping-induced 
enhancement of the PL yield and of its decay rate as seen in our intense excitation regime 
(see Fig. [3]), in which the doped charge density is negligible with respect to the photoinduced 
one, indicates that the radiative centers are probably not to be associated with bulk oxygen 
vacancies, Nb ions or other donor centers, in contrast to what often stated in the literature.— 
This is to be contrasted with the regime of low excitation intensity, in which doping-induced 
electrons are not negligible and doped samples do show a greatly enhanced yield and decay- 
rate of the PL.— i2&2£i2I The small, but significant, sample-to-sample fluctuations of PL yield 
and decay rate that we see and the presence of three different spectral components in the PL 
spectra having different thermal behavior concur to indicate that a significant role in the STO 
luminescence is played by unidentified intrinsic defects (other than standard bulk oxygen 
vacancies), most likely by providing the radiative centers responsible for the luminescence 
itself. A less likely alternative hypothesis is that the luminescence is a fully intrinsic process 
(e.g., taking place in STEs) and the defect role is that of introducing competing non-radiative 
decay channels which decrease the yield. In both cases, the responsible defects might be 
titanium interstitials,- 26 crystal dislocations,— or other defect complexes^ 1 ^ or, possibly, 
surface defects, as suggested by the strong PL enhancement seen in acid-etched samples 
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and in STO nanoparticles.— ^ On the possible role of the surface states, see also Ref. 
However, the PL spectral difference between Nb-doped and pure samples that we observed 
at high temperatures remains unexplained. 

Let us now move on to discussing the physical interpretation of the PL decay dynamics 
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that naturally emerges from our best model, i.e. C2PUB (vl or v2), with the help of Fig. [TT] 
which provides a schematic picture of the most important dynamical processes assumed to 
occur in the system. At time zero, the UV excitation will generate many electron-hole pairs 
in the CB and VB (process a), part of which will redistribute very quickly (within a time- 
scale of few picoseconds, negligible for our experimental resolution) in trapped states lying 
in the band-gap (horizontal solid lines in Fig. [TT]) . The initial conditions given in Table I 
are defined by this rapid redistribution process. Most likely the number and energy depth of 
these traps is not symmetrical between holes and electrons. Here, for definiteness, we assume 
that there are mainly trapped holes and free electrons, although the converse is equally 
possible. Trapped holes, or possibly trapped hole-electron pairs (i.e. trapped excitons) will 
form population 1 of model C2PUB while free electrons will correspond to population 2. 
It is also assumed that N\ <C N2, i.e., the trapped holes are only a small fraction of the 
total number, so that the free electrons and holes are approximately still balanced. These 
free pairs will then recombine directly in (defect-assisted or phonon-assisted) non-radiative 
bimolecular processes (process b in Fig. [TTT) controlled by rate constant 72. This provides 
the faster decay channel, entirely non-radiative and (approximately) independent of the N\ 
decay. The best-fit values of the time constant r 2 , combined with our estimate of the FDCF 
constant, leads to an estimated band-band recombination rate constant 72 — 10~ n s _1 cm 3 . 
This value is intermediate between the typical orders of magnitude found for indirect and 
for direct semiconductors,— which seems a plausible result. 

Thus far the two variants of model C2PUB are essentially identical. They are however 
distinguished by the interpretation attributed to population 1 and its two decay terms. 
The simplest (and hence more plausible) interpretation is that corresponding to variant 
C2PUBv2, described by the two processes under the label (c) in Fig. [TTJ In this case, the 
unimolecular term (with rate constant 1/ri) is simply taken to be a (non-radiative) thermal 
untrapping of holes (not balanced by trapping of free holes, except for very short times 
after excitation, because free holes decay much faster), while the coupling term proportional 
to N1N2 is taken to be a cross recombination between a trapped hole and a free electron, 
leading to the PL emission, with rate constant ji . Phonons will also be created in the 
emission (triggered by defect vibrational excitations, shown as dashed lines in Fig. [TTT) . 
thus explaining the PL Stokes shift and the spectral shape discussed in Sec. 3. The different 
spectral components (BL, GL and VL) are here ascribed to the existence of different kinds of 
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FIG. 11: (Color online) Schematic picture of the main electronic processes that may be postulated 
for interpreting the rate equations of our decay models C2PUB vl and v2: (a) - pump photon (UV) 
absorption and electron-hole pair generation; (b) - non-radiative direct electron-hole recombination; 
(c) radiative cross-recombination between a free electron and a trapped hole, with blue photon 
(BL) emission, and non-radiative thermal untrapping of holes; (d) - radiative crossed recombination 
between a free electron and a trapped exciton; (e) - radiative spontaneous annihilation of a trapped 
exciton. The constants 71,72, 1/ti are the dynamical rates of each process. 

trapping sites. We also note that in the intense excitation regime in which we have performed 
our measurements, it is likely that the initial number of trapped holes is saturated. This 
would explain quite naturally the initial conditions adopted in C2PUBv2. On the other 
hand, for lower excitation energies, a non saturated linear behavior Aq(0)~[7 should be 
resumed, thus explaining also the quadratic-to-linear yield crossover reported by Yasuda et 
al.M 

To interpret model C2PUBvl we must assume instead that population 1 corresponds to 
excitons (i.e. e-h pairs), rather than unpaired holes, that are either self-trapped or trapped 
near a structural defect. In this case, population 1 can decay either by spontaneous annihi- 
lation (process e in Fig. [IT]) or by "crossed" recombination of a hole of the trapped pair with 
a mobile electron, leading to the emission of one PL photon and the simultaneous freeing 
of the electron of the pair (process d in Fig. [TT]) . Both the unimolecular and the coupling 
terms would then be radiative with similar quantum yield, as assumed in model C2PUBvl. 
Although possible in principle, we believe that this C2PUBvl model scenario is somewhat 
less plausible than the former (the C2PUBv2 one), both because it is more complex and 
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because it is difficult to justify the very rapid creation of a population of trapped excitons 
as would be required by the initial conditions assumed in the model. 

In conclusion, the main results of the investigation reported in this Article are the fol- 
lowing: (i) we have presented strong evidence that the radiative centers involved in the 
blue luminescence cannot be associated with bulk oxygen vacancies or other donor impu- 
rities; (ii) we have shown, nevertheless, that a crucial role in the luminescence is played 
by other yet-unidentified intrinsic structural defects, such as dislocations, defect complexes, 
and possibly the surface; these defects likely provide the actual radiative centers; (iii) we 
have shown that the initial decay in the investigated excitation-intensity range is dominated 
by a bimolecular process, while trimolecular processes such as Auger are not significant; 
(iv) we have provided strong evidence that at least two separate interacting photoexcited 
charge populations are involved in the PL dynamics, which we interpret simply as mobile 
and defect-trapped charges. 
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